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ABSTRACT 



Context. During the course of a large spectroscopic survey of X-ray active late-type stars in the solar neighbourhood, we discovered 
four lithium-rich stars packed within just a few degrees on the sky. Although located in a sky area rich in CO molecular regions and 
dark clouds, the Cepheus-Cassiopeia complex, these very young stars are projected several degrees away from clouds in front of an 
area void of interstellar matter. As such, they are very good "isolated" T Tauri star candidates. 

Aims. We present optical observations of these stars conducted with 1-2 meter class telescopes. We acquired high-resolution optical 
spectra as well as photometric data allowing us to investigate in detail their nature and physical parameters with the aim of testing the 
"runaway" and "in-situ" formation scenarios. Their kinematical properties are also analyzed to investigate their possible connection 
to already known stellar kinematic groups. 

Methods. We use the cross-correlation technique and other tools developed by us to derive accurate radial and rotational velocities 
and perform an automatic spectral classification. The spectral subtraction technique is used to infer chromospheric activity level in 
the Ha line core and clean the spectra of photospheric lines before measuring the equivalent width of the lithium absorption line. 
Results. Both physical (lithium content, chromospheric, and coronal activities) and kinematical indicators show that all stars are very 
young, with ages probably in the range 10-30 Myr. In particular, the spectral energy distribution of TYC 4496-780-1 displays a strong 
near- and far-infrared excess, typical of T Tauri stars still surrounded by an accretion disc. They also share the same Galactic motion, 
proving that they form a homogeneous moving group of stars with the same origin. 

Conclusions. The most plausible explanation of how these "isolated" T Tauri stars formed is the "in-situ" model, although accurate 
distances are needed to clarify their connection with the Cepheus-Cassiopeia complex. The discovery of this loose association of 
"isolated" T Tauri stars can help to shed light on atypical formation processes of stars and planets in low-mass clouds. 

Key words, stars: fundamental parameters - stars: pre-main sequence - stars: formation - stars: kinematics - X-rays: stars 



1. Introduction 

Although most studies of very young low-mass stars have con- 
centrated on star-forming regions (SFRs) or young clusters, it 
has become increasingly evident that a considerable number of 
very young stars unrelated to prominent SFRs is present in the 
solar neighbourhood. Their existence far from interstellar ma- 
terial is difficult to reconcile with the standard picture of star 
formation. Different scenarios, such as "runaway" T Tauri stars 
(ISterzik & Durisen | 
clouds (Feigel sonl ll 1 



19951 fl998h or formation in small turbulent 
1996), have been proposed as possible expla- 
nations. Most of the comoving associations of nearly coeval stars 
(i.e., moving groups (MG)) discovered so far are located in the 
southern hemisphere and may be tied to the Sco rpi us-Centaurus 
assoc iation (see, e.g. JZuckerman & Song 2004 and lTorres et alJ 
2008). 



Single active stars in the field, selected on the basis of their 
high coronal emission, are mostly young stars with an age of 
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* Based on observations collected at the Observatoire de Haute 
Provence (OHP, France) and the Catania Astrophysical Observatory 
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** Figure [5] is only available in electronic form at 
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a few 100 Myr, i.e. in the zero-age main sequence (ZAMS), or 
even younger: post- T Tauri stars (PTTS) or T Tauri stars (TTS). 
iGuillout et af](l2009h conducted a large ground-based observing 
program designed to perform the phy sical characterization of 
* 1000 RasTyc stellar X-ray sources dGuillout et all 1 1999b by 
analysing their high-resolution optical spectra. Although the 
sky density of the youngest stars (identified thanks to their 
very high-lithium content) is more or less uniform dKlutschl 
2008), we discovered an unusual group of four lithium-rich 
stars towards the Cep heus-Cassiopeia (Cep-Cas) complex (see, 
e.g.. iKun etaD[2008l and references therein). However, their 
space distribution in a wide region devoid of dense CO clouds 
and their stellar properties make the relation between these stars 
and the Cep-Cas complex uncertain. Young stars in the sky 
region around the Cep-Cas complex are indeed not necessarily 
associated with these SFRs. A good example is represented 
by the nearby young visual binary HI P 115147 (V368 Cep ) 
and its reported comoving companion (Makar ov et al.l 12007). 
which are also projected towards the same region of the sky. 
HIP 1 1 5147 is currently classified as a very you ng "naked" 



PTTS dNations et alJ 119901: IChugainov etaiTll99lL 1 19931) . 20 
to 50 Myr old, located at 20 pc whose origin remains disputed 
although it is certainly unrelated to the Cep -Cas complex 
dMontes et alJl2001atlL6pez-Santiago et alJl2010b . 
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Table 1. Main data of our sources from the literature and dates of observations. 



Name 


RasTyc source a (2000) 


8 (2000) 


V T 


(B - V) T 


PSPC 


Li obs. 


Ho- obs. 




h m s 


o 1 II 






(ct/s) 






BD+78 853 


RasTycOOOO+7940 00 00 38.4 


+79 40 36.9 


10.35 


0.79 


0.10 


2001/11/10 


2001/11/20 


TYC 4496-780-1 


RasTyc0013+7702 00 13 40.5 


+77 02 10.9 


9.83 


0.70 


0.10 


2001/11/12 


2001/11/21 


TYC 4500-1478-1 


RasTyc0038+7903 00 38 06.0 


+79 03 28.5 


10.37 


1.04 


0.11 


2003/12/14 


2009/10/06 


BD+78 19 


RasTyc0039+7905 00 39 40. 1 


+79 05 30.8 


9.67 


0.77 


0.05 


2002/1 1/02 


2002/11/05 


Table 2. Johnson-Cousins and Stromgren photometric data. 


Name 


V R c 


h 


B-V 


U-B 


b-y 


iri[ 




BD+78 853 


10.29 (0.05) 9.77 a 


9.42" 0.62 (0.02) 


0.16(0.02) 


0.462 


0.164 


0.266 


TYC 4496-780-1 


9.78 (0.05) 9.22 a 


8.858 c 0.64 (0.03) 


0.15 (0.02) 


0.438 


0.211 


0.269 


TYC 4500-1478-1 


10.47 a (0.03) 9.96 d (0.03) 9.60 d (0.02) 0.87 (0.03) 


0.40 (0.02) 


0.572 


0.243 


0.400 


BD+78 19 


9.57 (0.03) 9.11 (0.03) 8. 


89 (0.02) 0.72 (0.03) 


0.09 (0.02) 


0.457 


0.218 


0.260 



Zach arias et al.l J2004I) b lMonet et al.l d2003h c iDroege et al J f200^) d Corrected from the contribution of the nearby companion. 



In the present paper, we concentrate on the four stars that we 
previously discovered and study their properties, evolutionary 
status, and kinematics. The observations and data reduction are 
outlined in Sect. [2] The determination of their stellar parameters 
as well as the level of chromospheric activity, lithium content, 
and age are briefly discussed in Sect. [3] We discuss the prop- 
erties of each star under scrutiny in Sect. |4] We conclude and 
summarize the results from the ongoing study of these sources 
in Sect.|5] 

2. Observations and data reduction 

Spectroscopic observations were conducted in 2001, 2002, and 
2003 (see Table [B at the OHP 1.52-m telescope using the 
Aurelie spectrograph to acquire spectra in both the Ha (AA « 
6500-6620 A) and lithium ( AA a 6650-6775 A) spectral regions 
at a resolution of R « 38 000. An echelle spectrum at R « 21 000 
of TYC 4500-1478-1 was also taken with the FRESCO spectro- 
graph of the 0.9 1-m telescope of the OAC in 2009. For details 
about the re duction of Aurelie and FRESCO spectra, the reader 
is referred to lGuillout et all (120091) . 

In 2008, complementary photometric observations were per- 
formed in the standard UBV system and Stromgren uvby filters 
with the 91 -cm telescope of the OAC and a photon-counting re- 
frigerated photometer equipped with an EMI 9893QA/350 pho- 
tomultiplier, cooled to — 15°C. For details a bout the reduction o f 
photometric data, the reader is referred to Fras ca et"ail (|2006). 
CCD images in Johnson-Cousins V, Rc, and Iq filters were 
also acquired for two stars using the OAC focal-reducer imag- 
ing camera. The reduction of these data was performed follow- 
ing standard methods and the VRqIc magnitudes were extracted 
from the corrected images by means of aperture photometry. 
Photometric data are summarized in Table [2] 

3. Astrophysical parameters, chromospheric 
activity, and lithium abundance 

High-resolution spectroscopic observations allow us to derive 
radial (RV) and projected rotational velocities (vsku), spec- 
tral type, luminosity class, and metallicity. They also enable us 
to identify spectroscopic binaries. Cross-correlation functions 
(CCFs) were computed to d erive radial velocit i es (Ta ble |4]l and 
we used the ROTFIT code dFrasca etal.ll2003l 120061) to evalu- 
ate effective the temperature (r e ff), gravity (logg), metallicity 



([Fe/H]), and v sin ; of our targets (TableO. We were also able to 
infer both the level of chromospheric activity (from the emission 
in the Ha line core) and the age (from the lithium abundance). 
The equivalent width of the lithium line (Wu) and the net equiv- 
alent width of the Ha line (WJSD were measured in the spectrum 
obtained after subtracting the non-active template b y integrat- 
ing the resi dual Ha and lithium p rofile. We refer to lFrasca et al.l 
(120061) and iGuillout et al.l d2009l) for a more detailed presenta- 
tion of the methods we used to derive lithium abundances, along 
with chromospheric and coronal luminosities (Table|3]l. 

As an independent estimate of the main astrophysical param- 
eters, we evaluated the effective temperature, radius, and metal- 
licity of these stars from our Stromgren photometry using the 
uvbyB algorithm dMoonll 19851) . which also provides an estimate 
of the color excess E(b - y). As seen in Table [3] the values of 
effective temperature derived from spectroscopy and Stromgren 
photometry agree with each other to within 200 K. The metallic- 
ity index 5mq is in substantial agreement with the spectroscopic 
value [Fe/H]«0 for all stars. Estimated radii and \ogg values 
excluded the possibility of lithium-rich giant stars. 

4. Discussion 

4.1. Star properties 

The high resolution spectra and CCFs of our stars are shown 
in Figs. Q]and |3E f° r both Ho- and Lii spectral regions. The 
characteristics of the stars, derived from the analysis of the 
spectra acquired in the Ha and lithium spectral regions, are now 
summarized below: 

BD+78853: Although both the Ha and lithium spectra 
display broad absorption lines, which are sometimes asymmet- 
ric, the automatic procedure failed to detect double peaks in 
any of the CCF The Ho- profile is seen in absorption, but the 
analysis found that the line is partly filled. The equivalent width 
of the lithium line is WLi = 0.18A, corresponding to a lithium 
abundance log N(Li) = 3.05, which implies that BD+78 853 is 
young. The rotational velocity derived from the line broadening 
decreases by a factor two or so from the Ho- (vsin i a 50 km s _1 ) 
to lithium spectra (v sin i « 28 km s _1 ). This could be indicative 
of a spectroscopic binary, but may also be the result of a different 



1 Available in electronic form only. 
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Table 3. Astrophysical parameters from spectroscopy and Stromgren photometry, lithium abundance, and chromospheric/coronal 
luminosities. 



Name 


1 eS 

(K) 


logg 1 


[Fe/H] 1 


1 eB 

(K) 


R 1 

(Rg) 


6m z 


E(b-yY 
(mag) 


WWlogJV(Li) 
(A)/ 


WSJ/login, 
(A) /(ergs" 1 ) 


logZ-x 
(ergs -1 ) 


BD+78 853 


5738 


4.19 


-0.03 


5750 


0.98 


0.032 


0.068 


0.18/3.05 


0.31/29.5 


30.6 


TYC 4496-780-1 


5844 


4.17 


-0.04 


5650 


0.99 


0.012 


0.024 


0.18/3.14 


3.90 / 30.6 


30.4 


TYC 4500-1478-1 


5160 


4.30 


-0.06 


5330 


0.88 


0.052 


0.099 


0.30/3.09 


0.97 / 29.6 


30.4 


BD+78 19 


5444 


4.14 


-0.11 


5630 


0.96 


0.003 


0.040 


0.22 / 3.00 


0.21/29.5 


30.2 



1 Derived from spectra 2 Derived from Stromgren photometry 



CCF shape in the two spectral domains. The radial velocities 
derived from Ha and Li i spectral regions are compatible within 
0.2 km s _1 at RV = -6.4km s . Our present data suggest that 
BD+78 853 is a single star that is rotating relatively quickly, but 
they do not exclude a spectroscopic binary system. More data 
are needed to settle this point. 

TYC 4496-780-1: The Ha line exhibits an absorption feature 
roughly centered on a strong asymmetric Ha emission profile. 
The Li 1/16707.8 line with W Li = 0.18A (log iV(Li) = 3.14) 
is stronger than the nearby calcium line, as in the case of 
BD+78 853. In both spectra, the lines are very broad and 
the code detected two Doppler-shifted dips (at -5.5 and 
45.2 km s ) in the Ha CCF with very different depths, sug- 
gesting a binary in which the two stars of the system have 
significantly different luminosities. In the lithium CCF, the 
footprint of a secondary component is clearly visible (although 
not detected by the automatic analysis) as an asymmetry 
on the red side of the main CCF dip, which is centered on 
about -4 km s _1 . We derived an RV ~25km s _I for the 
secondary dip, which would imply an RV variation of +2 and 
-20 km s _1 for the primary and secondary component, respec- 
tively, of this SB2 candidate. The projected rotational velocity 
of TYC 4496-780-1 is estimted to be v sin i = 37 +/ - 10 km s" 1 . 

TYC 4500-1478-1: The FRESCO spectrum displays an Ha 
profile markedly filled-in by emission (see Fig. 1). The lithium 
region shows very sharp absorption lines and exhibits a very 
strong Li i ,16707. 8 line with Wo =0.30 A corresponding to 
log/V(Li) = 3.09. The perfect fit of the CCF allows us to derive 
accurate values of radial velocity (RV = -8.8kms _1 ) and 
projected rotational velocity (vsinz s* 9kms~ 1 ). Unless we have 
observed a binary system at a conjunction, we consider TYC 
4500-1478-1 as a single star. 

BD+78 19: The Ha spectrum is typical of a mid-G type star 
without any remarkable characteristics. Both RV = -9.7km s _I 
and v sin i « 12 km s _1 were derived from the Gaussian fit of the 
CCF. Apart from the observed Li i /16707.8 deep absorption line 
(W Li = 0.22 A; log/V(Li) = 3.00), the lithium spectrum CCF is 
clearly not Gaussian and a fit of a rotational profile provides a 
closer match. The radial velocity determination from the lithium 
region is compatible with the former value but the rotational 
velocity is larger by a factor 2 or so. This suggests that, although 
the individual spectra appear as single-lined, BD+78 19 may be 
a binary system of mass ratio « 1 observed with Aurelie close to 
conjunction at both epochs. 

The chromospheric and coronal luminosities (see Table [3) 
were computed assuming that all stars are 15Myr old (see 
Table [4] for the adopted distances ofisMyr and Sect. 14.31 for its 



justification). We found X-ray luminosities Lx - 10 3a4 ergs~' 
for all sources (within 0.2 dex), which is typical of weak- 
line T Tauri stars (WTTS) in Taurus-Auriga-Perseus SFRs. For 
both BD+78 853 and BD+78 19, we found L Hn = 10 29 5 e rg s" 1 
a value slightly higher than that found by Guill out et alj ([2009) 
for their PTTS candidates and about an order of magnitude lower 
than that in old binary systems, in which the coupling of spin 
and orbital motions by tidal actions can maintain hig h chromo- 
spheric activity for a very long time (iFrasca et al.ll2006l) . The Ha 
luminosity of TYC 4496-780-1 is an order of magnitude higher 
than that of BD+78 853 and BD+78 19, but both the shape of the 
Ha profile and the IR excess (see Sect. 14. 2t infer that accretion 
is the primary cause of the observed emission. 

4.2. Spectral energy distributions 

The standard UBVRqIc photometry (Table com plemented 
with JHK S magnitudes from the 2MASS catalogue ( Cutri et all 
2003), allowed us to reconstruct the spectral energy distribution 
(SED) from the optical to the near infrared (IR) domain for all 
sources. 

We used the grid of NextGen low-resolutio n synthetic spec- 
tra, w ith log g = 4.0 and solar metallicity by lHauschildt et alJ 
(1999), to perform a fit to the SEDs. The effective temperature 
(Teff) was kept fixed to the value derived with ROTFIT (Table [3]) 
and we allowed the interstellar extinction (Ay), which affects 
the SED shape particularly at the shortest wavelengths, and 
the angular diameter, which scales the model surface flux 
to the stellar flux at Earth, free to vary. The Cardell Tet alJ 
(1989) extinction law with Ry - 3.1 was used. The optimal 
solution was found by minimizing the chi-square of the fit. 
The latter was performed only on UBVRIJ data, which are 
dominated by the photospheric flux of the star and are normally 
not appreciably affected by infrared excesses. The Ay values 
are in the range 0.06-0.2 mag, in agreement with the color 
excesses E(b - y), which were inferred from the Stromgren 
photometry. The rather low Ay values are typical of sources 
at « 80-200 pc that are not embedded in dusty nebulae, and 
consistent with their position far from dark clouds. The angular 
diameters (if), all on the order of 0. 1 mas, are reported in Table|4] 

Figure [2] displays the results of our fitting procedure. As 
is clearly evident, the SEDs are reproduced well by the syn- 
thetic spectrum blueward of the K s band, with the exception of 
TYC 4496-780-1, which displays a strong near- and far-IR ex- 
cess (IRAS fluxes). The lack of significant IR excess in three 
of our targets implies that these stars are likely to be either 
WTTS or PTTS. The slope of the SED between K s and 24yum, 
a =i —1.1, allows us to tentatively classify TYC 4496-780-1 as 
a Class II young infrared source, i.e. a T Tauri st ar surrounde d 
by an accretion disc, according to the definition of lLadaldl987l) . 



4 
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Fig. 1. Spectra (black line) in the Ha (upper panels) and lithium (lower panels) region of the Cepheus comoving stars, and the 
template spectra (red lines) of non-active lithium-poor reference stars broadened to the v sin i of the targets and Doppler-shifted 
according to their RV. 




2-6 3.8 4.0 4.2 4.4 3.5 4.0 4.5 5.0 5.5 6.0 3.6 3.B 4.0 4.3 4.4 3.6 3.B 4.0 4.2 4.4 

los(A [AD IoeCa [AH loa(A [All 1q E (A [All 



Fig. 2. Spectral energy distributions of the four stars (dots). The best-fit NextGen spectrum is overplotted with continuous lines in 
each box. 



Table 4. Parameters deduced from the SEDs and kinematics. The proper motions are from the TYCHO-2 catalogue. 



Name 


<P 


rfzAMS 




fi a cos 6 


Us 


RV 


l^ZAMS 


^ZAMS 


WzAMS 


V 15Myr 


VlSMyr 


Wl5Myr 




(mas) 


(pc) 


(pc) 


(mas yr" 


') 


(kms" 1 ) 




(km s" 1 ) 






(kms- 1 ) 




BD+78 853 


0.088 


105 


175 


23.5 


-0.4 


-6.5 


7.0 


-10.5 


-4.0 


13.9 


-13.9 


-6.0 


TYC 4496-780-1 


0.104 


100 


150 


22.4 


2.2 


-5.5 


6.5 


-9.9 


-1.9 


11.1 


-12.6 


-2.2 


TYC 4500-1478-1 


0.100 


75 


130 


25.1 


-3.4 


-8.8 


3.2 


-11.5 


-4.1 


8.9 


-14.6 


-5.3 


BD+78 19 


0.092 


90 


155 


23.2 


-2.8 


-9.7 


3.6 


-12.7 


-4.3 


9.8 


-16.3 


-5.5 



This classification is consistent with its Ho- emission profile. 
Similar Ha line profiles are often observed in accreting TTS 
( App enzeller et al] |2005). According to the concept of the mag- 
netospheric accretion model, the large line width (several hun- 
dred km s~') is indicative of large-scale gas flows, while the in- 
verse P Cygni profile traces the gas infall (iMuzerolle et alj2 001 : 
iKurosawa etal1l2006l) . 



4.3. Kinematics and relation to Cep-Cas complex. 

The spatial distribution of BD+78 853, TYC 4496-780-1, TYC 
4500-1478-1 and, BD+78 19 is show n in Fig. [3] overplotted on 
the extinction map (Av) pub lished by [ Dobashi et al. (20051), m 
a similar way to Fig. 2 of iTachihara et alJ d2005l) . The posi- 
tion of the nearby young visual b inary HIP 1 15147 and its co- 
moving companion, discovered by Makarov et al.1 (120071) . is also 
plotted together with the "isolated" young stars discovered by 
ITachihara et al.l d2005l) . The closest prominent star-forming re- 
gion is the Cep-Cas complex, which surrounds our stars. They 



are located in particular, at Galactic latitudes 14° < b < 17°, 
on the west side of the Cepheus flare, a region of signifi cant ex- 
tinctio n outside the Galactic plane first recognized by [Hubble 
(Il934h . Across the Cepheus flare region. lKur3(ll998h determined 
the distances to several dark clouds and concluded that the inter- 
stellar matter is concentrated at 200, 300, and 450 pc (although 
the whole Cep-Cas complex extends to 800 pc or s o) separated 
rough ly at regular intervals along Galactic latitudes. Olan o et al.l 
(2006) suggested that the Cepheus flare forms part of a large 
shell, 50 pc in radius, ex panding at 4 km s' 1 , that encloses an 
old supernovae remnant dGrenier et al.lll989h . whose center can 
be located at / = 120°, b = 17°. The Ce pheus flare I s a reg ion of 
active star formation as demonstrated bv lKun et alJ (|2009), who 
discovered about 70 new PMS stars in the star-forming clouds. 

All the four stars investigated by ourselves are projected in 
front of the void located between the dark clouds LDN 1333 
(west), LDN 1241-1251 (east), and north of LDN 1259-1262. 
A group of seven stars, including TYC 4500-1478-1, located 
in this void and displaying typical T Tauri characteristics, was 
identified by Tach ihara et al.l (120051) . For a few stars, including 
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Fig. 3. Spatial distribution of our "isolated" comoving TTS (big 
diamonds) overplotted on a map of dark clouds enhancing in- 
terstellar material. The "isolated" young stars previously discov- 
ered in the CO void region are shown with square symbols, while 
the young "naked" PTTS HIP 115147 and its comoving com- 
panion are indicated by small asterisks. The major clouds in the 
vicinity are also labelled. 



TYC 4500-1478-1, they found a lithium equivalent width close 
to that often measured at the Pleiades upper envelope, implying 
a PTTS or ZAMS evolutionary status for them. Although they 
consider these stars to be associated with the Cep-Cas complex, 
only on the basis of their proximity to this region, and adopt a 
distance of 200 pc for them, their relation to the cloud complex 
remains uncertain. 

The distance of our stars is crucial in determining their 
origin and a possible association with the Cep-Cas complex. 
Unfortunately, the Tycho parallaxes are useless and one must 
rely on photometric distance estimates. To cover all the possi- 
bilities, we estimated for each star two distances that we as- 
sume to be a lower and upper limit. The lower limit was cal- 
culated for the hypothesis that the star is on the ZAMS. The 
star radius of a ZAMS star with the same was adopted and 
the distance was deduced from the angular diameter <p reported 
in Table [4] The upper limit was estimated assuming a mini- 
mum age and the corresponding radius was computed from the 
ISiess et ail d2000h evolutionary tracks. An age of 5 Myr would 
place stars at distances mostly exceeding 250 pc, i.e. farther 
than the closest boundary of the Cepheus flare shell in that 
direction. This would contradict our photometric observations 
that no stars experience significant interstellar extinction. An 
age older than w 20-30 Myr would contradict the observational 
findi ng that, at this age, stars display no significant NIR ex- 
cess (iHillenbrand & FEPS Spitzer Legacy Science Team.l l2006: 
Mey er et alj 120061) . A good compromise is 15 Myr that we 
adopted in the following. The lower and upper limits to the dis- 
tances (c/zams and disMyr, respectively) are reported in Table|4] 

As seen in this table, the ZAMS assumption would place 
our stars roughly within lOOpc of the Sun ruling out their 
association with the Cep-Cas complex whose CO clouds are 
farther away by nearly a factor of two. Needless to say that 
this assumption is also hardly compatible with the detection of 
the IR excess from an accretion disc around TYC 4496-780-1. 
On the other hand, the 15 Myr assumption would place them 
roughly at the distance of LDN 1333, LDN 1261, and LDN 1228 




-30 -20 -10 10 

V„ (km s") 



Fig. 4. U- V kinematic diagram for our comoving stars (labelled 
as in Table |5). Blue circles and red triangles refer to the ZAMS 
and 75 Myr assumptions, respectively. The average velocity 
components (dots) of some young SKG and those of some late- 
type stars members of these young SKG are also plotted (square, 
triangle, circle, upside -down triangle, and U symbols for the 
IC 2391 supercluster, Pleiades, Castor, UMa moving groups, and 
Hyades supercluster, respectively). The loci of the young-disc 
(YD) and old-disc (OD) populations are also marked. The posi- 
tion of HIP 1 15147 is plotted as a filled square. 

(see Fig.5 of iKun et alJ 120081) and can thus partly reconcile 
their possible association with at least the nearest clouds of the 
Cepheus Flare. 

Additional crucial information to elucidate their origin ten- 
tatively can be derived from kinematics. Because errors in the 
distances are the major uncertainties (proper motions and radial 
velocity errors are smaller than or equal to 2 mas and 1.5 km s _1 
respectively), we computed the heliocentric space velocity com- 
ponents, UVW, in a left-handed coordinate system, both with 
distances corresponding to ZAMS radii (I/zams VzamsWzams) 
and adopting the rfisMyr distances (J/i5MyrVi5MyrWi5Myr)> i- e - as ~ 
suming that all stars are located in proximity to the closest clouds 
of the Cep-Cas complex. The RV values measured o n our spec- 
tra and Tycho-2 proper motions (lH0g et alJ l2000t) have been 
used. The positions of our four stars and HIP 115147 in the U- 
V kinematical diagram are shown in Fig. [4] together with the 
mean pos ition of the major stell ar kinematics groups (SKG) dis- 
cussed in Mon tes et al.l d2001bl) . namely the IC 2391 superclus- 
ter (~ 50 Myr), the Pleiades SKG (~ 100 Myr), the Castor SKG 
(~ 200 Myr), the Ursa Major (UMa) group (~ 300 Myr), and the 
Hyades supercluster (~ 600 Myr). Figure |4] readily shows that, 
whatever the distance is, all our stars share the same kinematics 
(within a few km s _I ) proving that they form a homogeneous 
MG with the same origin. We also note that both distances and 
kinematics seem to rule out an origin in common origin with 
HIP 115147. 

To assess their membership to the aforementioned SKGs on 
more obje ctive gro u nds, w e used the probabilistic approach de- 
scribed in Klutsch (2008), i.e. we computed the probability P 
that a star with heliocentric velocities U, V, and W has a Galactic 
motion compatible with a given SKG. Independently of the dis- 
tance assumption, the quantitative analysis excludes any possi- 
ble association with the oldest UMa and Hyades MGs (P « 0), 
as expected from the high-lithium photospheric abundance of 
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Table 5. Membership probability P (%) of our four stars to ma- 
jor MGs computed with the ZAMS (left) and 15Myr assumption 
(right). 



# 


Name 


IC 2391 


Pleiades 


Castor 


1 


BD+78 853 


0.0/4.7 


15.9/36.9 


59.6/41.7 


2 


TYC 4496-780-1 


0.0 / 0.3 


8.9 / 22.7 


33.4/28.2 


3 


TYC 4500-1478-1 


0.0/0.1 


11.3/50.5 


33.2/32.8 


4 


BD+78 19 


0.0/0.3 


17.4/65.9 


28.6/20.0 



our stars. Table suggests a possible link with the Castor MG, 
which is, however, again not supported by the very high-lithium 
abundance measurements, hardly compatible with a 200 Myr old 
MG. With a less than 5% probability, the IC 2391 MG ori- 
gin is highly improbable. We are left with the Pleiades SKG 
for which the probabilities rise from « 15 % for the ZAMS as- 
sumption to 40 % for th e 15 M yr assumption. However, as re- 
called by M akarov et al.1 (120071) . this stream includes isolated 
stars, groups, and associations of diverse ages in the range from 
1 Myr to about 200 Myr, which by itself, does not shed light on 
the origin of our stars. 

We finally consider a possible link with the Cep-Cas 
complex. Projected on the sky, the closest dark clouds to our 
stars are LDN 1251 and LDN 1 241 at both 5-8° east of them 
and about 300pc from the Sun dKun et al.ll2008l) . The average 
proper motions of LDN 1251, as estimat ed thanks to its cloud 
members discovered by iKun et all d2009l) . are p a = 1.9+ 1 .7 
and fig = -0.5 + 1.3masyr _1 . The angular distance of 5-8° 
corresponds to 25-40 pc at the cloud distance. With a proper 
motion of about 21 masyr 1 relative to LDN 1251, nearly 
aligned in the east-west direction and corresponding to a 
velocity of about 30km s at the cloud distance, they should 
have travelled for 0.8-1.3 Myr from their parental cloud. The 
closest cloud on the approachin g side of t he expanding shell 
(namely LDN 1228, d 180pc. Kmienill2008l) is about 10° 
east of our stars (« 30 pc at the cloud distance). By completing 
the same calculations for LDN 1251/1241, we found an escape 
velocity (v esc ) from LDN 1228 of about 22 kms~' and a travel 
time of about 1.3 Myr. These estimates do not allow us to 
exclude that these stars are "runaway" objects that originated in 
LDN 1251, LDN 1241, or LDN 1228. Nevertheless, although a 
close encounter wi th massive stars c a n res ult in high velocity 
"runaway" stars, ISterzik & Durisenl dl998l) found that 60% 
of the runaways have speeds larger than 3kms~' but that 
this fraction steeply decreases, becoming smaller than 5 % 
for v esc >10kms _1 . Both the high escape velocites and the 
common space motions cast doubt on the "runaway" hypothesis 
being applicable. 

The more plausible explanation of how these "isolated " 
TTS formed is given by the "in-situ" model dFeigelsonl 19961). 
This s cenario would support the conclusion of Tachihara et al. 
(12001 although the lack of accurate distances does not al- 
low us to draw firm conclusions about the origin of this group 
of stars. The advent of GAIA mission, with its unprecedented 
astrometric precision, will certainly shed light on this open 
issue. Similar conc lusions have been rea ched in the Taurus 
(Mag azzu et all 1997k Chamaeleo n-Musca dMizuno et al.l 1998) 
and Lupus (Tachih ara et al.l [2001) SFRs although no accreting 
TTS have been found outside these SFR's cores. 



5. Conclusion 

We have reported the discovery of four comoving very young 
stars located in a region devoid of dark matter and molecular 
clouds. On the basis of its Ha profile and infrared excess, one of 
the stars investigated has been tentatively classified as a Class II 
young infrared source, i.e. a TTS surrounded by an accretion 
disc. Owing to the distance uncertainty, we cannot assert that 
these stars are completely unrelated to the Cep-Cas star-forming 
region. However, their kinematics prove that they form a homo- 
geneous comoving group of stars with the same origin. Their 
off-cloud positions imply that they are very good "isolated" TTS 
candidates. The "runaway" hypothesis is highly improbable be- 
cause of their kinematical properties. Our conclusions raise a 
question about the applicability of the "in-situ" star-formation 
scenario to very low-mass cloud environments. TYC 4496-780- 
1 could be a TW Hya analog, although a little bit older, and its 
comoving companions may represent the peak of the iceberg of 
a young loose association. 

Acknowledgements. We are grateful to the OHP night assistant staff in con- 
ducting our Key Program, and those of the OAC observatories for their sup- 
port and help with the observations. We thank Luigia Santagati of INAF-Catania 
for the English revision of the text. This research made use of SIMBAD and 
VIZIER databases, operated at the CDS, Strasbourg, France. This publica- 
tion uses ROSAT data. Part of this work was supported by the Universidad 
Complutense de Madrid (UCM), the Spanish MICINN, Ministerio de Ciencia 
e Innovation under grant AyA2008-00695 and the Comunidad Autonoma de 
Madrid, under PRICIT project S-2009/ESP-1496 (AstroMadrid). 

References 

Appenzeller, I., Bertout, C, & Stahl, O. 2005, A&A, 434, 1005 
Cardelli, J. A., Clayton, G. C, & Mathis, J. S. 1989, ApJ, 345, 245 
Chugainov, P. E, Lovkaya, M. N., & Petrov, P. P. 1991, Information Bulletin on 

Variable Stars, 3623, 1 
Chugainov, P. E, Petrov, P. P., & Lovkaya, M. N. 1993, Izv. Krym. 

Astrofizicheskoj Obs., 88, 29 
Cutri, R. M., Skrutskie, M. E, van Dyk, S., et al. 2003, The IRSA 2MASS All- 
Sky Point Source Catalog, NASA/IPAC Infrared Science 
Dobashi, K., Uehara, H., Kandori, R., et al. 2005, PASJ, 57, 1 
Droege, T. E, Richmond, M. W., Sallman, M. P., & Creager, R. P. 2006, PASP, 

118,1666 
Feigelson, E. D. 1996, ApJ, 468, 306 

Frasca, A., Alcala, J. M., Covino, E., et al. 2003, A&A, 405, 149 
Frasca, A., Guillout, P., Marilli, E., et al. 2006, A&A, 454, 301 
Grenier, I. A., Lebrun, E, Arnaud, M., Dame, T. M., & Thaddeus, P. 1989, ApJ, 
347, 231 

Guillout, P., Klutsch, A., Frasca, A., et al. 2009, A&A, 504, 829 

Guillout, P., Schmitt, J. H. M. M., Egret, D., et al. 1999, A&A, 351, 1003 

Hauschildt, P. Ft., Allard, E, & Baron, E. 1999, ApJ, 512, 377 

Hillenbrand, L. & FEPS Spitzer Legacy Science Team. 2006, Bulletin of the 

American Astronomical Society, 38, 910 
H0g, E., Fabricius, C, Makarov, V. V., et al. 2000, A&A, 355, L27 
Hubble, E. 1934, ApJ, 79, 8 

Klutsch, A. 2008, Ph. D Thesis, ULP and Strasbourg Observatory 
Kun, M. 1998, ApJS, 115, 59 

Kun, M., Balog, Z., Kenyon, S. J., Mamajek, E. E., & Gutermuth, R. A. 2009, 
ApJS, 185, 451 

Kun, M., Kiss, Z. T, & Balog, Z. 2008, Handbook of Star-Forming Regions, ed 

B. Reipurth (San Francisco: ASP), Vol. 1, 134 
Kurosawa, R„ Harries, T. J., & Symington, N. H. 2006, MNRAS, 370, 580 
Lada, C. J. 1987, IAU Symposium, eds. M. Peimbert & J. Jugaku, 1 15, 1 
Lopez-Santiago, J., Montes, D., Galvez-Ortiz, M. C, et al. 2010, A&A, 514, 

A97 

Magazzu, A., Martin, E. L., Sterzik, M. E, et al. 1997, A&AS, 124, 449 
Makarov, V. V., Zacharias, N., Hennessy, G. S., Harris, H. C, & Monet, A. K. B. 

2007, ApJ, 668, L155 
Meyer, M. R., Hillenbrand, L. A., Backman, D., et al. 2006, PASP, 118, 1690 
Mizuno, A., Hayakawa, T, Yamaguchi, N., et al. 1998, ApJ, 507, L83 
Monet, D. G., Levine, S. E., Canzian, B., et al. 2003, AJ, 125, 984 
Montes, D., Lopez-Santiago, J., Fernandez-Figueroa, M. J., & Galvez, M. C. 

2001a, A&A, 379, 976 
Montes, D., Lopez-Santiago, J., Galvez, M. C, et al. 2001b, MNRAS, 328, 45 



P. Guillout et al.: Discovery of "isolated" comoving T Tauri stars in Cepheus (RN) 



1 



Moon, T. 1985, Ap&SS, 117, 261 

Muzerolle, J., Calvet, N., & Hartmann, L. 2001, ApJ, 550, 944 

Nations, H. L., Heckert, P. A., Mooney, C. B., et al. 1990, BAAS, 22, 1254 

Olano, C. A., Meschin, P. I., & Niemela, V. S. 2006, MNRAS, 369, 867 

Siess, L., Dufour, E., & Forestini, M. 2000, A&A, 358, 593 

Sterzik, M. F. & Durisen, R. H. 1995, A&A, 304, L9 

Sterzik, M. F. & Durisen, R. H. 1998, A&A, 339, 95 

Tachihara, K., Neuhauser, R., Kun, M., & Fukui, Y. 2005, A&A, 437, 919 

Tachihara, K., Toyoda, S., Onishi, T., et al. 2001, PASJ, 53, 1081 

Torres, C. A. O., Quast, G. R., Melo, C. H. F., & Sterzik, M. F. 2008, Handbook 

of Star Forming Regions, Volume II, 757 
Zacharias, N., Monet, D. G., Levine, S. E., et al. 2004, A&AS, 36, 1418 
Zuckerman, B. & Song, I. 2004, ARA&A, 42, 685 



P. Guillout et al.: Discovery of "isolated" comoving T Tauri stars in Cepheus (RN), Online Material p 1 





-200 -100 100 200 -200 -100 100 200 -200 -100 100 200 -200 -100 100 200 

Radial velocity (km/s) Radial velocity (km/s) Radial velocity (km/s) Radial velocity (km/s) 



Fig. 5. Cross-correlation functions (cross symbols) of the Cepheus moving groups TTS candidates BD+78 853 
(RasTyc 0000+7940), TYC 4496-780-1 (RasTyc 0013+7702), TYC 4500-1478-1 (RasTyc 0038+7903), and BD+78 19 
(RasTyc 0039+7905) (from left to right) computed from their Ha (upper panels) and lithium (lower panels) spectra. The 
continuous red line displayed in each panel shows the best fit, either with a Gaussian or a rotational profile. The residuals are 
displayed with a thin black line in the bottom of each panel. 



